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Abstract 



The existence of a two-partial Rayleigh wave coupled to an electrical 
P] , field in 2nim piezoelectric crystals is known but has rarely been investigated 

O ' analytically. It turns out that the Z-cut X-propagation problem can be fully 

solved, up to the derivation of the secular equation as a polynomial in the 
squared wave speed. For the metallized (unmetallized) boundary condition, 
^ I the polynomial is of degree 10 (48). The relevant root is readily identified 

OO ' and the full description of the mechanical and electrical fields follows. The 

^H ■ results are illustrated in the case of the superstrong piezoelectric crystal, 

\^ [ Potassium niobate, for which the effective piezoelectric coupling coefficient 

—il ' is calculated to be about 0.1 
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1 Introduction 

This article prolongs and complements papers by the present authors [Ij and 
by others [21 |3l HJ [5l El [TJ |8] where the propagation of a Shear- Horizontal (SH) 
surface acoustic wave, decoupled from a two-partial Rayleigh surface acoustic 
wave, was considered for piezoelectric crystals. Those papers examined situ- 
ations (cuts, propagation directions) where the interaction between acoustic 
fields and piezoelectric fields concerns the SH wave exclusively and not the 
Rayleigh wave, which remains purely elastic. In the present paper, the situ- 
ation is reversed: the interaction occurs solely between the electric field and 
the mechanical displacement lying in the sagittal plane (the plane containing 
the direction of propagation and the direction of attenuation), leading to a 
piezoacoustic two-partial (elliptically polarized) Rayleigh surface wave. 

The properties of a two-partial Rayleigh surface wave complement those 
of a SH surface wave and one wave's loss is the other's gain. Hence SH surface 
waves are particularly suited for immersed crystals (liquid sensing, biosensors, 
etc.) because the mechanical displacement is polarized horizontally with 
respect to the interface, which leads to low loss of acoustic power in the 
fluid; conversely, two-partial Rayleigh surface waves are used extensively for 
non-destructive surface evaluation [9] and for free surface sensors [TO], because 
their propagation is highly sensitive to anything present on the interface 
which might perturb their vertical displacement. To take but one example 
it is possible, using Rayleigh surface waves, to design a mass microbalance 
with a mass resolution of 3 picograms[llj. 

This context reveals the importance of studying the analytical properties 
of such waves. The cuts allowing for the propagation of two-partial Rayleigh 
waves coupled to an electric field were identified and classified by Maer- 
feld and Lardat [12] ; these waves were also investigated numerically [I3l [H] 
and experimentally [15], as is best recalled in the textbook by Royer and 
Dieulesaint [16] (see also Mozhaev and Weihnacht [8] for pointers to more 
recent contributions.) In general, the problem treatment however falls short 
of a full analytical resolution, and the wave speed is usually found from a 
trial- and-error procedure which goes back and forth between the propaga- 
tion condition and the boundary condition, until a certain determinant is 
minimized to a required degree of accuracy [I9] (alternatively, Abbudi and 
Barnett [H] proposed a numerical scheme based on the surface-impedance 
matrix.) The present paper shows that a secular equation can be derived 
explicitly as a polynomial of which the wave speed is a root, for the Z-cut 



X-propagation problem. 

This feat is achieved by use of some fundamental equations (II) satisfied 
by the 6-vector whose components are the mechanical displacements and 
tractions and the electrical potential and induction at the interface. Albeit 
powerful, the method based on the fundamental equations has one drawback 
because the polynomial secular equation possesses several spurious roots. 
Hence for the metallized boundary condition (III. A), it is a polynomial of 
degree 10 in the squared wave speed, and for the unmetallized boundary 
condition (III.B), it is a polynomial of degree 48! Nevertheless, finding the 
numerical roots of a polynomial is almost an instantaneous process for a 
computer. Also, it is expected that among all the 10 or 48 possible roots, 
one gives exactly the surface wave speed. Consequently, that root satisfies 
the boundary condition exactly, whereas none of the spurious roots does. 
Once the relevant root is thus properly identified, all the quantities of in- 
terest follow naturally: the attenuation coefficients, the depth profiles, the 
electromagnetic coupling coefficient, etc. Here, the method is applied to the 
superstrong piezoelectric crystal. Potassium niobate KNO3, for which the ef- 
fective electromagnetic coupling coefficient for the piezoacoustic surface wave 
is found to be about 0.1. 



2 Basic equations 



2.1 Constitutive equations and equations of motion 

Consider a piezoelectric crystal with two mirror planes (orthorhombic 2mm, 
tetragonal 4mm, or hexagonal 6mm). For this type of crystal, the elasto- 
piezo-dielectric matrix[20] is of the form. 
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• • • • 

• • • • 



(1) 



Now consider the Z-cut, X-propagation of a surface acoustic wave that 
is, a motion with speed v and wave number k where the displacement field 
u and the electric potential (p are of the form, 

{u,(P}{x,,X2,X3,t) = {U{kx3),^{kx3)}e''^^^-'"\ (2) 

(say), with 

C/(oo) = 0, (^(oo) = 0. (3) 

Here the xi, X2, x^ axes are aligned with the crystallographic axes, and the 
crystal occupies the X3 > region. 

It follows from the constitutive equation Eq. ([T]) that the tractions aij 
and the electric induction Di are of a similar form, 

{aij, Di}{xi, X2, X3, t) = tk{tij{kx3), di{kx3)}e'''^'''~'"\ (4) 

(say) with ^22 = -ic23U-^ - 1632^' + C12U1, d2 = -ie24U2, 

til = -icisf/g - iesiifi' + ciiUi, 

tl3 = -ic^bU'i + C55f/3 + ei5V3, 

te = -^casf/s - iessv^' + Ci3?7i, 

t23 = —10441/2 + C4QU2, ti2 = —iC4QU2 + CqqU2, 

di = -iei^U'i + eiiUi + Cigf/g - env?, 

4 = -^633^^3 + ^e33V5' + 631 f/i, (5) 

where the prime denotes differentiation with respect to kx^. Also, the surface 
wave vanishes away from the interface, so that 

ty(oo) = 0, di{oo) = 0. (6) 

The classical equations of piezoacoustics, aijj = pui^tt, Di^i = (where p 
is the mass density of the crystal), reduce to 

- hi + iti3 = -pv'^Ui, - ti2 + it'23 = -pv'^U2, 

- ti3 + it^jg = -pv^Us, -di + id'^ = 0. (7) 

Clearly, the second equation Eq. ([7]) 2 involves only the function U2 and is 
decoupled from the three others, which involve the functions Ui, f/3, and ip. 
It reads: C44U2 — {cqq — pv^)U2 = 0. A simple analysis shows that there are no 



functions U2 solution to this second-order ordinary differential equation such 
that f/2(oo) = and ^23(0) = — 2C44f/2(0) = 0, except the trivial one. Hence, 
the piezoelastic equations, coupled with free surface boundary condition, lead 
to plane strain: U2 = 0, which in turn leads to (generalized) plane stress: 
ti2 = t23 = d2 = by Eqs. I^^as- 

Now the remaining constitutive equations and piezoacoustic equations 
can be arranged as a first-order linear differential system. It develops as: 



^' = tNt 



(8) 



where (using the notation of Ting [27]) 
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(9) 



Lothe and Barnett[21j established the explicit expressions for the compo- 
nents of the real matrix N in the general case (general anisotropy, general 
piezoelectricity). In the present context, they are given by 



X = pv\ 6 

■56 = 615/055, 

nee = I/C55, ^^22 = ^33/6'^ n24 = €33/6 

rj = Cii- [ci3(e33Ci3 + 2631633) - 033631]/^^ 



= £33633 + 633, 

^4 = (£33613 + 631633)/^^ r2 

2 „ „ 7^2 



= (613633 - 631633)/^^ 

n44 = -633/'^^, 

Ai = -(en -^615/655). (10) 



2.2 General solution 

The solution to the linear system with constant coefficients Eq. ([H]) is of 
exponential form. Indeed, taking ^ as: C,e^^'^^^ where (^ is a constant vector 
and q a decay coefficient, leads to the eigenvalue problem: {N — ql^jC, = 
where Iq is the 6x6 identity matrix. Hence, g is a root (with positive 
imaginary part, to ensure decay) to the propagation condition: det{N — 
qIq) = 0, which is a cubic for g^. 



0J4q'^ + uj2q^ 



LOq = 0, 



(11) 



where 

uji = n22X + nQQ{X -rj) - n^^jj, + 2r2SQ + 2r4, 

U2 = {{X - C55)[e33(X - Cii) - 631] 

- eii[X(c33 + C55) - C11C33 + Ci3(ci3 + 2C55)] 

- ei5[X(ei5 + 2e3i) + 2633) - 2633011 + 2(ei5 + 631)013)]} 

/[c55(e33C33 + e33)], 
Uo = -{X - Cii)(Xeii - 6^5 - eiiC55)/[c55(e33C33 + 633)]. (12) 

Here of course, it must be realized that the propagation condition Eq. ( II ip 
can be solved for q only once the speed of the surface wave (and hence 
X = pv"^) is known. The next subsection and the next section show how 
X can be found as a root of the secular equation. Once X is known, the 
propagation condition gives six roots, out of which only three are kept: gi, 
g2, Q'3 say, the three roots with positive imaginary roots ensuring exponential 
decay (if for a given X, the propagation condition fails to deliver three such 
roots, then no surface wave can propagate at speed a/a/p.) 

Let C"*-, C^, C,^ be the corresponding eigenvectors: NCJ" = qiC (^ = 1; 2, 3), 
obtained for example as the third column of the matrix adjoint to N — qil^. 
Explicitly they are of the form, 

(^* = [aj,6j, — Ci,C55fi,C55gi,eQ — hi] , (13) 

£33 £33 

where the non-dimensional quantities a^, g^, hi contain only even powers of 



q, and the non-dimensional quantities Cj, fi, gi contain only odd powers of q: 

es^csso-i = -(esscss + 633)?^ + [e33(X - C55) - £11033 - 2ei5e33]g^ 

+ eii(X-C55) -ei5, 
esscss&i = [(ci3 + C55)e33 + 633(615 + e3i)]g^ + [(ci3 + C55)eii + 615(615 + 631)]^, 

ei5C33Q = -[633(615 + 631) + 633(613 + 655)]^^ + [6i5(X + 613) + 631 (X - 655)]^, 

633633655/1 = -655(633633 + 633)?^ 

+ [e33655(-^ + 613) - 633(611655 + 615 + 615631) + 633(615613 + 631655)]^^ 
+ [(611655 + 615) (X + 613) + X6i563i]g, 

6336336555'* = -655(633633 + 633)?^ + [(633613 + 633631) (X - 655) 

+ 6i5633(X - 613) + 633(611655 + 615 + 6i563i)]g^ 
+ 6i3[eii(X-655) -615], 
61560633/^1 = 615(633633 + 633)5^ 

- [e336l5(^ + 613) - 633(611655 + 615 - 615631) + 611(631633 - 633613)]^^ 
+ 63l[eii(X-655)-6?5]. (14) 

Then the general solution to the equations of motion Eqs. ([H]) is 

^(kxs) = 7i(^ie^'='?i^^ + 72C^6^'='?2^« + 73C=^6^'^'?3^^ (15) 

where 71, 72, 73, are constants. 

Depending on the type of boundary conditions, a given homogeneous 
system of three linear equations for 71, 72, 73 is derived. The corresponding 
determinantal equation is the boundary condition. In general for surface 
waves, the interface X3 = remains free of tractions: ^31(0) = ^33(0) = 0. 
From these two equations, 72 and 73 can be expressed in terms of 71 as 

72 ^ /35'i - fl93 73 ^ /i5'2 - /2gl .^g^ 

7i f293 - /35'2 ' 7i f293 - hg2 ' 

To sum up: first the speed of the surface wave must be computed as a 
root of the secular equation (Section III), obtained thanks to the fundamen- 
tal equations presented below (Section II. C). Next the appropriate decay 
coefficients are computed as roots with positive imaginary parts from the 
propagation condition Eq. fITT]) . Then it must be checked that the boundary 
condition (Section III) is indeed satisfied. If it is, then the complete solution 
is given by Eqs. (0), (EDi, (USD, (USD- 



2.3 Fundamental equations 

Now some fundamental equations are presented, from which the secular equa- 
tion is found. Their derivation is short and is given in Refs. [221 |23l [2l]; they 
represent a generalization to interface waves of works by Currie[25] and by 
Taziev|2S] for elastic surface waves (see also Ting [27] for a review.) They 
read 

'0 J3I 



^(0) • M("^|(0) = 0, where M^"^ 







AT", 



(17) 



and n is any positive or negative integer. By computing the integer powers 
N^ of N [aX n = —2, —1, 1, 2, 3, say), it is a simple matter to check that the 
6x6 matrix M^'^^ is symmetric and that its form depends on the parity of 
n. Hence M^'^'> is of the forms. 
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(18) 



when n = —2, 2, and when n = —1, 1, 3, respectively. 



3 Z-cut, X-propagation 

3.1 Metallized boundary condition 

For metallized (short-circuit) boundary conditions, the mechanically free in- 
terface 0:3 = is covered with a thin metallic film, grounded to potential 
zero, and so 
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(19) 



Two possibilities arise for the roots with positive imaginary part of the 
bicubic Eq. ( ITT]) . Either (a) qi = iqi {qi > 0) or (b) gi = — g^, ^3 = iq^ 
(gs > 0). In case (a), it is clear from Eqs. f fT4|) that at, gi, hi, are real 
numbers and that bi, Ci, fi are pure imaginary numbers. Then, separating 
the real part from the imaginary part in the third, fourth, and fifth lines in 
Eq. f lTI?]) 9. it is found that [71,72,73]^ is parallel to a real vector. It follows 
from Eq. f lTI?]) i that ^(0) is of the form 



C(0) = f/i(0)[l,^a2,0,0,0,/3i 



(20) 



where ia2 := f^3(0)/f/i(0) is pure imaginary (02 is real) and (3i := d3{0)/Ui{0) 
is real. In case (b) a slightly lengthier study shows that ^(0) is also of 
this form (see Ting|27] and Destrade|21] for proofs in different, but easily 
transposed, contexts.) 

Now substituting this expression Eq. ( 120|) for ^(0) into the fundamental 
equations Eq. (TT7|l i leads to a trivial identity when n = —2,2, and to the 
following set of three equations when n = — 1, 1, 3, 
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Note that the components of the 3x3 matrix and of the right hand-side 
column vector above are easily computed from their definition Eq. flT7|) 9: for 



instance, M, 
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Cramer's rule applied to the system above reveals that 2/3i = A2/A, 
Pf = A3 /A (where A is the determinant of the 3x3 matrix in Eq. (12 Ij) 
and A2, A3 are the determinants of the matrix obtained from this matrix by 
replacing the 2"'^ and 3'''^ columns by the right hand-side column in Eq. (J2T]) . 
respectively) and so, that 



A^ - 4AA3 



0. 



(22) 



This is the explicit secular equation for the speed of a two-partial Rayleigh 
piezoacoustic surface wave propagating in a metallized mm2 (or 4'm''m', or 
6mm) crystal. 

Its expression is too lengthy to reproduce here but has been obtained 
using Maple. It turns out that the secular equation is a polynomial of degree 
10 in X. Note also that the solution to the system Eq. fl^Tl) for the unknown 



a| plays no role in the final expression of the secular equation. Hence the 
equation is only valid in the presence of piezoelectric coupling through the 
solutions of Eq. fl2T]) for 2/3i(= 2d3{0)/Ui{0)) and for (3f, and it does not cover 
the Rayleigh cubic function for purely elastic surface waves in orthorhombic 
crystals. Moreover, in the present context the in-plane piezoacoustic surface 
wave is entirely decoupled from its anti-plane counterpart (which does not 
exist, as seen in II. A) , and so the secular equation Eq. fl2^ cannot cover 
the case of a Bleustein-Gulyaev SH wave. Hence in many respects, this new 
secular equation is unique and stands alone, with no link whatsoever with 
previously established secular equations. 

Selecting the correct root or roots out of the 10 possible given by the 
secular equation is quite a simple matter. First the root X must be real 
and positive; then it must be such that the propagation condition Eq. (1 11 1) 
written at X yields three roots gi, q2, q^ with positive imaginary part; finally 
it must be such that the boundary conditions Eqs. flT9|) 9 are satisfied, that is 



1 



{Qi - (l2){q2 - q3){q3 - Qa) 



Cl C2 C3 

/i /2 /a 
9i 92 93 



0. (23) 



For Potassium niobate[28] (KNbOa, 2mm), the relevant constants are the 
following. Elastic constants (10^^ N m~^): cn = 2.26, C13 = 0.68, C33 = 1.86, 
C55 = 0.25; Piezoelectric constants (C m^^): 615 = 5.16, 631 = 2.46, 633 = 4.4; 
Dielectric constants (10"^^ F m"^): en = 34eo, €33 = 24eo, eo = 8.85416 ; 
Mass density (kg m~^): p = 4630. The secular equation has six complex roots 
and four real positive roots in X; out of these four, only two are such that 
the propagation condition yields suitable attenuation coefficients; out of these 
two, only one is such that the boundary condition is verified. The numerical 
values for the wave speed Vm (say) and for the attenuation coefficients are 
listed on the second line of Table I. A 8 digit precision is given, although the 
calculations were conducted with a 30 digit precision; the left hand-side in 
Eq. fl23|) was found to be smaller than 5 x 10^^^. The complete solution is 
found by taking the real part of the right hand-side in Eq. (|2]) and Eq. (^. 
Specifically, the mechanical displacement Ui is in phase quadrature with the 
mechanical displacement U3 and the electric potential, 

ui = Uii^kxs) cos/i;(xi — vt), 

U3 = Usi^kxs) sin k{xi —vt), 

(p = 0{kxs) smk{xi ~ vt), (24) 



10 



where Ui := 5R{t/i}, f/3 := Q^jt/s}, ip := Q{4>} are the amplitude functions. 
Figure 1 shows their variations with the scaled depth X3/A, where A = 27r//c is 
the wavelength. The vertical scaling is such that f/i(0) = lA. The axes of the 
polarization ellipse are along the xi and X3 axes. At the interface, f/i(0) > 0, 
f/3(0) < 0, and |f/3(0)| > |f/i(0)|, so that the major axis of the ellipse is 
along X3 and the minor axis is along xi; there, the ellipse is spanned in the 
retrograde sense with time. The ellipse becomes more and more oblong with 
depth, and is linearly polarized at a depth of about 0.174A. Further down 
the substrate, it becomes elliptically polarized again, but is now spanned in 
the direct sense. It is circularly polarized at a depth of about 1.183 A, and 
then again linearly polarized at a depth of about 0.987A. 

3.2 Unmetallized boundary condition 

For the unmetallized (free) boundary condition, the free surface is in contact 
with the vacuum (permeability: eo), and so [I], 
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Similarly to the previous case, the form of ^(0) can be found, whatever 
the form of the qi is. Namely, 



^(0) = <^(0)[^a2,/3l,l,0,0,^eo] 



(26) 



where ia2 = Ui{0)/ip{0) is pure imaginary (^2 is real) and (3i = f/s(0)/(y9(0) 
is real. Substitution into the fundamental equations Eqs. flT7|) at n = —2, 2 
leads to the trivial identity. At n = —1, 1, 3 it leads to 
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(27) 



which are three equations for two unknowns 0:2 and /3i. Formally, solving 
two equations and substituting the result into the third equation yields the 
secular equation. Note however that these equations Eqs. f l7r|) are nonlinear 
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(quadratic) in the unknowns. Their resolution is somewhat lengthy, although 
possible as is now seen. 

First take advantage of the identity M23 = N^e = to solve Eqs. fl27|) at 
n = 1 for Pf: 



r(l) 



r(i) 



r(i). 



= [/U - n^^el + 2r2eoa2 + (^ - X)al)]/X. 



Next, solve Eqs. f l27j) at n = — 1, 3 for /3i 



2A = [Mg) + e^M^e") + 2eoM£)a2 + M^;^)^^' + m£)/3?]/M 



,(n) 
23 5 



(28) 



(29) 



n = —1,3. Now square both sides and substitute the expression for f3l just 
obtained to derive two polynomials of fourth degree in 02. Having 0:2 as 
a common root, these two polynomials have a resultant equal to zero, a 
condition which is the explicit secular equation for the speed of a two-partial 
Rayleigh piezoacoustic surface wave propagating in an unmetallized 2mm (or 
4mm, or 6mm) crystal. 

Of course, the resulting polynomial is rather formidable, here of degree 
48 in X according to Maple. Nevertheless, finding numerically the roots of 
a polynomial is a quasi- instantaneous task for a computer. For instance in 
the case of KNbOa, it is found that there are 10 positive real roots in X to 
the polynomial, out of which 6 yield three attenuation factors with positive 
imaginary part. Out of these 6, only one satisfies the boundary conditions 
Eqs. ([25]) 2, that is 



1 



(gi -?2)(g2-g3)(g3-g4) 



/l /2 /s 

9i 92 93 

hi - ici /i2 - ic2 hs - ic3 



(30) 



Hence, at the speed Ki (say) and attenuation factors listed on the third line 
of Table I (obtained with a 40 digits precision), the determinant in Eq. (1301) 
was found to be smaller than 5 x 10~^^. Note by comparison of the second 
line and the third line of Table I that when the free surface is metallized, the 
wave propagates at a slower speed, and is slightly more localized, than when 
the surface is unmetallized. Figure 2 shows the variations of the amplitude 
functions with the scaled depth X3/X in the unmetallized boundary conditons 
case. The depth curves are similar to those in the metallized case, with the 
differences that the boundary condition forces the electrical potential to be 
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about 0.596 V at the interface, and that the nature of the polarization elhpse 
changes at depths which are shghtly less than the corresponding depths with 
metallized boundary conditions. 

Finally, recall that it is usual to take the quantity 2(Vu — Kn)/Vu as a 
measure of the crystal's ability to transform an electric signal into an elas- 
tic surface wave by means of interdigital electrode transducers although, as 
proved by Royer and Dieulesaint [12] , the demonstration is far from obvious. 
This quantity is often referred to as the effective piezoacoustic coupling coeffi- 
cient for surface waves and is expected to be positive. In the present example 
of KNbOs, the speeds of the second column in Table I give a value of 0.1037, 
far greater than the corresponding values [TB] for GaAS, Bii2Ge02o, ZnO, 
and CdS, and more than twice that [22] for LiNbOa. Note that Mozhaev and 
Weihnacht [17] reported negative values for this quantity corresponding to 
special cuts and propagation direction in KNbOa. 
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Table I. Wave speed (m s ^) and attenuation coefficients for a two-partial 
piezoacoustic surface wave in KNbOa. 

metallized 3762.50953 ± 0.39191249 + i 0.49991830 i 3.10691826 

unmetallized 3968.28624 ± 0.39840475 + i 0.45628905 i 3.07008806 
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Figure 1: Depth profiles of the mechanical displacements (A) and the elec- 
tric potential (V) for the piezoelectric Rayleigh wave in KNbOs, Z-cut X- 
propagation with metallized boundarj^conditions. 




Figure 2: Depth profiles of the mechanical displacements (A) and the elec- 
tric potential (V) for the piezoelectric Rayleigh wave in KNbOs, Z-cut X- 
propagation with unmetallized boun(|^y conditions. 



